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Abstract 

Classically, percolation critical exponents are linked to the power laws that characterize 
percolation cluster fractal properties. It is found here that the gradient percolation power 
laws are conserved even for extreme gradient values for which the frontier of the infinite 
cluster is no more fractal. In particular the exponent 7/4 which was recently demonstrated 
to be the exact value for the dimension of the so-called " hull" or external perimeter of the 
incipient percolation cluster keeps its value in describing the width and length of gradient 
percolation frontiers whatever the gradient value. Its origin is then not to be found in the 
thermodynamic limit. The comparison between numerical results and the exact results 
that can be obtained analytically for extreme values of the gradient suggests that there 
exist a unique power law from size 1 to infinity which describes the gradient percolation 
frontier, this law becoming a scaling law in the large system limit. 

Spreading of objects in space with a gradient of probability is most common. From 
chemical composition gradients to the distribution of plants which depends of their solar 
exposure, probability gradients exist in many inhomogeneous systems. In fact inhomogeneity 
is a rule in nature whereas most of the systems that physicists are studying are homogeneous 
as they are thought to be more simple to understand. In particular phase transitions or critical 
phenomena are studied in that framework, the simplest being percolation transition^. 

In this work, the opposite situation, a strongly inhomogeneous system is studied. We 
report the discovery that the validity of some percolation exponents can be extended to 
situations which are very far from the large homogeneous system limit. This is found in the 
frame of gradient percolation, a situation first introduced in the study of diffusion front[||, ||. 
Surprisingly these exponents, up to now believed to pertain to large systems, are also verified 
in a limit that could be called the small system limit (S.S.L.) where some of the gradient 
percolation properties can be computed analytically. This letter presents a report concerning 
the 2D square lattice. The same results and discussion have been obtained for the triangular 
lattice and will be published elsewhere, together with the details of the analytical calculations 
for both lattices. 

The gradient percolation (G.P.) situation is depicted in FIG. |[ The figure gives an 
example of a random distribution of points on a lattice with a linear gradient of concentration 
in the vertical direction. It is a 2D square lattice of size L g x L, where each point (x, y) is 
occupied with probability p(x) = 1 — x/L g (x being the vertical direction in the figure). In 
gradient percolation there is always an infinite cluster of occupied sites as there is a region 
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where p is larger than the standard percolation (S.P.) threshold p c . There is also an infinite 
cluster of empty sites as there is a region where p is smaller than p c . The object of interest is 
the G.P. front, the external limit (or frontier) of the infinite occupied cluster. It is constituted 
by the sites which belong to the occupied cluster and are first-nearest neighbours with empty 
sites belonging to the infinite empty cluster. It is shown in grey in FIG. |l[ This front is a 
random object with an average position xj, a statistical width a f and a total length Nf. In 
so far that the G.P. front and the S.P. external perimeter (often called hull) have the same 
geometry their fractal dimension was first conjectured to be exactly equal to 7/4 in Q. This 
result was then demonstrated heuristically by Saleur and DuplantierQ and very recently it 
was proved mathematically by Smirnovpj. 

The early G.P. studies were focussed at finding its relation with standard percolation. Let 
first recall the definitions. For < x < L g , nt(x) is the mean number of points of the front 
lying on the line x per unit horizontal length. It measures the front density at distance x. 
The length Nf, the position Xf and the width Of of the front are then defined in terms of the 
n f (x) by 
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It was found that the mean front was situated at a distance where the density of occupation 
was very close to p c or p(xf) ~ p c . This was verified numerically with such precision that the 
gradient percolation method is now often used to compute percolation thresholds ||, Q, ||, ||, 10]. 
It was also found that: 

1. The width cry depends on L g through a power law Of oc (Lg) v ^ 1+U ^ where v = 4/3 is 
the correlation length exponent jjj in dimension d = 2 so that c/ = (L 9 ) 4 / 7 . The width 
Of was also shown to be a percolation correlation length. 

2. Secondly it was found that the front was fractal with a dimension Df, numerically 
determined, close to 1.75. The front length followed a power law Nf oc {L g ) aN with 
a N = {D f - l)u/(l + u). 

3. But also, it was numerically observed that the sum of these two exponents was very 
close to 1. If true this meaned that v/(l + u) + {Df — 1)^/(1 + v) = 1 or Df = 1 + 1/v 
= 7/4. This is how it was conjectured in Q that Df = 7/4. 

In that sense the G.P. power laws were thought to be linked to the S.P. exponent v and 
to the fractality of the percolation cluster hull. Up to now, these facts were considered to be 
strictly valid only in the large system limit. 

However, if true, and we know now that 7/4 is the exact value, there follows an intriguing 

Df i 

relation, namely a j is exactly proportional to Lg. This means that the number of surface 
particles within the correlation length is exactly proportional to L g . This is particularly 
striking for diffusion fronts. Diffusion of particles from a source results in a concentration 
gradient and an associated G.P. situation. In that frame, the above result means that, if 
L g particles have diffused on a vertical row, there is on average the same (or a constant 
fraction of) number of particles on the correlated surface (surface content of a box with a 
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Figure 1: Gradient percolation front. Particles are distributed at random with probability 
p(x) = 1 — x/Lg. The occupied sites are in black and the last line of connected occupied sites 
is the gradient percolation front shown in light gray. This situation corresponds also to a 
diffusion situation where the front is called the diffusion front. Top: L = L g = 500. Bottom: 
L = 50 and L g = 12. The black window has an horizontal width equal to <r/. It contains 
approximately L g points. 
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lateral size equal to the statistical width). This fact seems a priori to have nothing to do 
with scaling, percolation and the thermodynamic limit. From this point of view it is possibly 
the consequence of a conservation law and if such a conservation exists, it should apply also 
for extreme gradients corresponding to L g of a few units. In particular it should apply to 
the very extreme L g = 1, 2 and 3 for which exact values of Xf, Nf and af can be calculated 
analytically. 

For L g = 2 or L g = 3, given a site on a line x, one can describe all the configurations 
such that the point belongs to the front, and compute their probability. For example, for 
L g = 2 all the occupied sites on the line x = 1 belong to the front, and a site on the line 
x = will belong to the front if at least one of its three neighbours on the line x = 1 is 
empty. Thus we get in this case, nj(l) = 1/2 and re/(0) = 1 — (1/2) 3 = 7/8. With the same 
kind of arguments, we can make the computations for L g = 3 but the geometry of connected 
sets is more complex as there are more configurations to consider. This however can be done 
exactly. One obtains: for L g = 1, aj = 0; 
• for L g = 2: 

11 4 , 2^7 



N f /L = —, x f = — and a f 
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The problem is then to compare the numerical G.P. laws to these exact values. As will 
be shown, the numerical results verify the above power laws with such precision that the 
question arises of the existence of a simple mathematical power law extending from L g = 1 
to infinity. To try to answer this question we proceed in 2 steps. First we test these laws on 
the numerical results obtained for L g between 4 and 50 for the square and triangular lattices 
by searching the best numerical power laws followed by the width. Considering arbitrary 
exponents a between 1.6 and 1.9, we study a" as a function of L g between 4 to 50. For each 
a value, there is a best line aj = a a {L g + b a ) fitting the numerical aj. The introduction 
of the term b a is justified by the fact that when a power law is verified for large systems 
it includes always the possibility that a small (as compared to the system size) term could 
contribute but in a negligible manner. But here the size itself is small or very small. On 
the other hand, one should remark that for L g = 1 the width is strictly so that some 
negative value of b a should be present. In the next step, the mean error d(a), defined by 
d(a) 2 = (1/47) X^ig=4(° 'f(L g ) a — a a (L g + b a )) 2 , is measured numerically as a function of a. 
The results are shown in FIG. |2[ There is a clear minimum for a = 1.75, showing that this 
exponent gives the best power law fit. Once the best fit with the empirical data is made one 
has the best values for the parameters a and b: a = 0.297 and b = —1.09. Note that b should 
be strictly equal to —1 in order to obtain a null width for the trivial case L g = 1. 

An other verification of the extreme G.P. power laws can be obtained from the study of 
the front length or of the quantity (Nf/L) 7 / 3 as a function of L g . In FIG. ||, the diamonds 
represent the values of (Nf/L) 7 / 3 and the best linear fit has equation Y = c(L g + d) with 
c = 0.845 and d = 0.88. This shows indeed that the exponents 4/7 and 3/7 can be used down 
to the steepest gradients for which the frontier is no more fractal. 

One can then extrapolate the af values to the case L g = 1, 2 and 3. The results are 
given in Table I. One observes that the numerical extrapolations correspond to the exact 
values with an apparent good precision. However no firm conclusion can be drawn without 
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Figure 2: Determination of the best exponent value for the square lattice: a = 1.75. The 
same result has been obtained for the triangular lattice. 
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Figure 3: Numerical results for the square lattice. The circles resp. diamonds represent 
respectively a 7 / 4 and (Nf/L) 7 / 3 . 
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discussion of the numerical uncertainties. The values shown on FIG. |3| are averaged over 50 
trials on a length L = 10 5 . As the fit occurs through a power law it is difficult to give a 
confidence interval for the coefficients a and b (obtained from a least square linear regression 
on the values of c^ 4 )- 

In order to obtain a better control on the numerical precision of a and b we made extensive 
computations of the two cases L g = 4 and L g = 5 with 100 trials on a length L = 5.10 5 . Doing 
so we obtain the mean values with their standard deviation: <jy(4) 7 / 4 = 0.8658 db 0.0009 and 
for <7f(5) 7//4 = 1.1610 ± 0.0013. Thus if we compute the equation of the line a(L g + b) which 
interpolates the two points (4, cr/(4) 7 / 4 ) and (5, cr/(5) 7/4 ), we obtain a = 0.2952 ± 0.0022 and 
b = —1.066 ± 0.041. This last result shows that the value —1 is not incompatible with b and 
its statistical error. Given the numerical values for L g = 4 and 5, we can also get extrapolated 
values for af for L g smaller, together with their confidence interval. The result (see Table I) is 
that the predicted values are very close to the exact ones. For (Nf/L) 7 / 3 , in the same way we 
obtain a linear interpolate of the values for L g = 4 and 5, with coefficients c = 0, 893 ± 0.014 
and d = 0.427 ± 0.074. 

At this point one could conclude that the numerical results are compatible with the 
existence of a single mathematical law for the width dependence, this law working from 
L g = 1 to infinity. Note however that the quality of the random number generator should 
intervene. As long as a mathematical proof has not been given, it is thus impossible for us to 
conclude on the exact values of the coefficients a and b. 

The above discussion bears on the validity of the numerical results. But the question of a 
unique mathematical power law can also be studied from the point of view of the exact results 
only. One can note that, remarkably, the values of the term b in the fit of the width are close 
to —1. If the above power law exists from L g = 1 to infinity it suggests that the real value of 
b is exactly -1 as the width is null in the trivial case L g = 1. As we have exact values, one can 
compute the equation of the line y = a(L g + b) defined by the two points (2, oy(2) 7//4 ) and 
(3,cr/(3) 7 / 4 ). One obtains a = 0.294 and b = -1.055. These values are close to the values 
obtained from the above numerical fit (a = 0.297 and b = —1.09), and here again b is close 
but not equal to —1. 

Why is there a small mismatch with the simplest law? The answer to this question is two- 
fold. First, it is possible that the exact power law is not valid for L g = 1 or both L g = 1 and 2. 
These cases could be "anormal" as for these values there is no Grossman- Aharony effect [11]. 
It is then possible that these two cases do not enter the unique mathematical law. Secondly 
the small discrepancy could be related to the fact that the frontier definition considers only 
the occupied sites. It gives to these sites a privilege role whereas one should also consider 
the frontier of the empty cluster. In fact this is not new in percolation studies || |8| where it 
was shown that the barycenter between the frontier of the occupied cluster and the frontier 
of the empty cluster was a more natural object. It notably permitted better computations 
of the percolation threshold. We have studied the statistical width of the local barycenter 
which can also be computed exactly for L g = 2 or 3. The results show the same behavior as 
described above i.e. a b value close but not equal to —1. The question remains then open to 
define the nature of the geometrical object which really would show a corresponding b value 
equal to —1. The answer has certainly interesting consequences for the S.P. problem itself. 

The above results have a simple practical consequence. Given an irregular non-fractal 
interface, for instance the grey line shown in Fig. 1, one can determine if it belongs to 
gradient percolation by measuring its statistical width a and the average value of Nf/L. 
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L g = 2 


L g = 3 


exact cjj 





0.4810 


0.7266 


07 (4-50) 


-0.13 


0.47 


0.72 


o s (4-5) 


0.106 


0.478 


0.726 


8a f (4-5) 


0.023 


0.005 


0.002 


exact Nf/L 


1 


1.3750 


1.6120 


N f /L (4-50) 


1.24 


1.48 


1.68 


N f /L (4-5) 


1.109 


1.393 


1.615 


5{N f /L) (4-5) 


0.017 


0.009 


0.004 



Table 1: Comparison between exact and extrapolated results. The data (4 — 50) resp. (4 — 5) 
correspond to extrapolated values from the respective ranges (4 — 50) resp. (4 — 5) (see text). 
5 is the confidence interval. 



Through the results described above one can find L g from the measured a f. Then one can 
check if Nf/L satisfies approximately the G.P. power law as a function of L g . This gives an 
intrinsic way to check if a given interface is of the gradient percolation type without previous 
knowledge of the gradient. This is important as there exist cases where irregular fronts like 
corrosion fronts belong to G. P. even though the gradient which built the interface is no more 
present and only the interface remains [12, 13]. 

In summary, it has been shown that the classical power laws of gradient percolation can be 
extended to extreme gradients with the same fractal exponents although the systems present 
no fractal geometry. Four comments can be drawn on these results. First, extreme gradient 
situations can be found in diffused contacts between materials. The fact that the contact 
geometry can be described by the same set of exponent whatever the width of the diffused layer 
is certainly of help to understand the properties of these contacts. There is an intrinsic method 
to find if a given rough interface belongs to gradient percolation without knowledge of the 
gradient. Secondly our results imply that there exists a conservation law which stipulates that 
the length of the correlated frontier is strictly proportional to the gradient length. Thirdly, 
the fact that the same exponent has been found for the square and the triangular lattice even 
for extreme gradients, in other words for small systems, suggest that universality here is not 
related to the neglect of the microscopic details of the interactions. Here there is no coarse 
scale and still universality is verified. Finally, the fact that the exponents 4/7 and 3/7 are valid 
down the smallest L g values (or the steepest gradients) suggests that these exponents play 
the same type of role here that the exponent 1/2 intervening in the fluctuations of the sum of 
independent identical random variables. In that last case the exponent apply to any number 
of random variables starting from 1, 2 or 3 up to infinity. The exponents 4/7 and 3/7 may 
then play a stronger role that critical exponents which only exist only in the thermodynamic 
limit. 

The Centre de Mathematiques et de leurs Applications and the Laboratoire de Physique 
de la Matiere Condensee are "Unite Mixte de Recherches du Centre National de la Recherche 
Scientifique" no. 8536 and 7643. 
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